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NOMENCLATURE 


Rry^Hah  Letter  Symbols 

A  surface  area 

a  exponent 

Cp  specific  heat  at  constant  pressure 
D  diameter 

F  Fahrenheit 

a  mass  flow  Telocity,  p  u 

0 

g  conrerslon  factor,  4.17  x  10 

o 

H  enthalpy 

h  convective  heat  transfer  coefficient 

J  conversion  factor,  778 

k  thermal  conductivity 

L  tost  section  length 

P  absolute  pressure 

p  gage  pressure 

q  heat  transfer  rate 

q'  heat  transfer  rate  per  unit  langth 

q"  heat  flux 

q*»*  heat  generation  per  unit  volume 
R  Ranklne 

tR  recovery  factor 

r  Radial  coordinate 

T  tea5)erature,  °R 

t  temperature, 

u  bulk  velocity  of  gas 

w  flow  rate 

X  axial  coordinate 

Subscripts  (refers  to) 

a  axial 

av  average 

aw  adiabatic  wall 

b  evaluated  at  bulk  temperatiire 


sq  ft 

Btu/(lb)(®F) 

ft 

lb/(sq  ft) (hr) 

(lb)(ft)/(lbp(hr^) 

Btu/lb 

Btu/(hr)(sq  ft)(®F) 
(ft)(lbj)/Btu 
Btu/(hr)(ft)(°F) 
ft 

Ibj/sq  ft 
Ibj/sq  ft 
Btu/hr 
Btu/(hr)(ft) 
Btu/(br)(sq  ft) 
Btu/(hr)(cu  ft) 


ft 


ft/hr 

Ib/hr 

ft 


e 

entrance  to  test  section 

f 

evaluated  at  film  temperatvire , 

Tj,  =  0.5(VTp«^ 

o 

stagnation  conditions 

V«) 

outside  s'srface  of  test  section  tube 

r 

radial 

T 

total 

t 

test  section  tube 

w 

inside  surface  of  test  section  tube 

0.4 

evaluated  at  reference  temperature, 

NOTE 

!  The  lack  of  a  subscript  on  gas  properties 

indicates  bulk  properties 

Greek  Letter  Symbols 

heat  flux  parameter 

q"  Ip  ■^!/c  V  T 

^  r  w  V  p  w  w 

a(  )  change  of 

P 

viscosity 

lb/ (hr) (ft) 

P 

density 

Ib/cu  ft 

T 

shear  stress 

Ib^/sq  ft 

P(  ) 

function  of 

Iton«dimensional  Parameters 

®Nu 

Nusselt  number 

hD/k 

«Rr 

Prandtl  nianber 

pC  /k 

P , 

«Re 

Reynolds  ntmber 

puD/p 

T 

w 

Tb 

temperature  parameter 

z 

D 

axial  position  parameter 

f 

friction  factor 

2 

2 

£H_ 


2^0 
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ii  research  Drogram.  entitled  "Hoat  Transfer  to  a  Gas  at  ttj.gh  'I'exaperature" 
is  being  conducted,  at  the  Stanford  Universitv  Nuclear  Technology  Laboratorj?- 
under  contract  inth  the  United  States  Atomic  Energy  Commissionj,  and  this 
report  presents  the  results  of  the  second  phase  of  that  program.  The  results 
of  the  first  phase  are  given  in  Technical  Report  No,  24’7-l  (2)  , 

Local  convective  heat  transfer  coefficients  and  local  and  average 
friction  factors  >jere  measured  for  nitrogen  floi-jing  through  an  electrically 
heated,  smooth,  round  tube.  Data  are  reported  for  a  range  of  Reynolds 
numbers  from  1800  to  90,000,  surface  temperatures  to  1^55°F,  gas  temperatures 
to  780°F,  and  heat  fluxes  to  138,000  Btu/(hr)(sq  ft). 

It  was  found  that  the  local  heat  transfer  data  could  be  correlated  hy 
the  equation, 

for  Reynolds  numbers  above  10,000, 

There  was  considerable  variation  in  local  friction  factor  along  the 
test  section;  however,  the  values  determined  near  the  exit  (but  far  enough 
from  the  exit  to  eliminate  end  effects)  were  correlated  by  the  relation, 


=  2  log 


for  Reynolds  numbers  above  10,000  with  fluid  properties  evaluated  at  the 
average  of  the  wall  and  the  gas  bulk  temperature.  Average  friction  factor 
data  were  also  correlated  by  this  equation  for  Re3molds  numbers  above 
20,000  but  with  fluid  properties  evaluated  at  bulk  temperatvire ,  At  lower 
Reynolds  numbers,  the  experimental  friction  factors  fell  considerably  above 
these  predicted  by  the  equation. 


Ntitnber  in  parentheses  refers  to  the  number  of  the  reference  appearxng  in 
the  List  of  References. 


2.  INTRODUCTION 


2.1  DEFINITION  OF  HEAT  TRANSFER  COEFFICIENTS 

In  the  analysis  of  convective  heat  transfer  processes  between  surfaces 
and  fluids,  it  is  common  practice  to  employ  a  convective  heat  transfer 
coefficient,  h,  to  relate  the  heat  flux  between  the  surface  and  the  fluid, 
q^,  at  a  point,  to  the  difference  between  the  surface  temperature,  T^, 
and  a  characteristic  temperature  of  the  fluid  at  that  point.  For  incom¬ 
pressible  pipe  flow  ,  the  local  heat  transfer  coefficient  is  ordinarily  based 
on  the  mean  stream  temperature.  For  high-speed  compressible  flow,  however, 
it  is  more  appropriate  to  base  h  on  the  adiabatic  wall  temperature,  "^aw’ 
\diich  is  the  temperature  that  an  insulated  tube  wall  would  assume  at  the 
same  R^olds  number  and  Prandtl  number.  Vfi.th  any  other  method  of  defining 
h,  it  would  be  possible  for  h  to  have  both  negative  and  positive  values. 
Moreover,  it  has  been  shown  (15)  that  the  value  of  h  defined  in  terms  of 
(T  -  T  )  is  the  only  one  leading  to  an  experimental  correlation  which  is 
independent  of  the  rate  of  heat  transfer.  Accordingly,  the  local  heat 
transfer  coefficient  is  defined  by  the  equation 

q«  =  h  (T  -  T  )  1 

^  '■  w  aw"^* 


The  adiabatic  wall  temperature  is  commonly  expressed  in  terms  of  the 
recovery  factor 


T  =  T  +'<R  (T  -  T)  2 

aw  'o' 

Shapiro  (15)  indicates  that,  to  a  good  approximation, 

^  =  N^^  3 

for  laminar  flow  and 

for  turbulent  flow.  Thus  for  air  the  recovery  factors  for  the  laminar  and 
turtulent  boundary  layers  are  approximately  0.85  and  0.90,  respectively. 


-  2  - 


The  use  of  a  local  heat  transfer  coefficient  is  perfectly  clear.  However, 
in  those  cases  where  primary  interest  is  focused  on  the  total  heat  transfer 
rate,  and  detailed  tempera tvire  and  heat  flux  distributions  are  of  no  particular 
concern,  it  has  been  fotuid  convenient  to  define  average  heat  transfer 
coefficients  (in  contrast  to  local  values  as  in  equation  1)  by 


av''  w,av 


aw,av' 


la 


In  the  case  of  heat  transfer  processes  for  fluids  flowing  through  circular 
tubes,  the  heat  flux  and  temperature  averages  are  properly  calculated  by 
dividing  the  integrated  values  along  the  heated,  or  cooled,  section  by  the 
total  length  of  the  section. 


2,2  DEFINITIONS  OF  FRICTION  FACTORS 

Pressure  loss  for  internal  flow  of  fluids  through  channels  is  most  con¬ 
veniently  related  to  the  flow  parameters  through  the  use  of  a  friction  factor. 
The  local  friction  factor,  f,  may  be  defined  as 


T 


Momentum  considerations  for  the  steady  flow  of  fluids  through  circular 
tubes  lead  to  the  follovdng  expression  for  the  apparent  wall  shear  stress, 

T^,  in  terms  of  quantities  which  may  be  determined  experimentally: 

T  5p,Ei!  I  5 

This  expression,  vdiich  neglects  the  change  in  elevation,  assmes  that  the 
static  pressure  is  uniform  across  the  flow  cross-section  and  neglegts  the 
effect  of  changing  velocity  profile  shape  on  the  momentum  flux,  •  The 

pressure  is  measured  by  means  of  taps  on  the  test  section,  while  the  accelera¬ 
tion  term  is  computed  using  the  bulk  gas  temperature  and  the  known  flpw  rate* 
An  average  friction  factor  may  be  defined  as 


f 


av 


T 

w,av 


7 
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or  equivalently  as 


f 


av 


8 


The  average  apparent  wall  shear  stress,  over  a  distance  from  station 

1  to  station  2  may  be  calculated  from  the  equation, 


It  is  often  useful  for  correlation  pxirposes  to  define  a  "film"  friction 
factor: 


10 


where  is  the  fluid  density  evaluated  at  some  reference  tenperature  as 
later  discussed  in  Section  3* 


-  - 


3.  RELEVANT  CORRELATIONS  AND  TECHNIQUES 


A  number  of  empirical  and  semi-analytical  methods  for  predicting  heat 
transfer  coefficients  and  friction  factors  have  arisen  from  investigations 
of  forced  convection  heating  of  gases  flowing  turbulently  through  circular 
tubes,  and  a  brief  review  of  the  more  common  methods  is  given  below. 
Particular  attention  is  focused  on  those  correlations  which  are  pertinent 
to  the  case  of  high  tenq^eratures  and  large  temperattire  differences  between 
the  surface  and  gas. 

3.1  THERMAL  ENTRANCE  REGION* 

Investigations  (3, ^,^,7, 10,14, 23)  of  convective  heat  transfer  have 
shown  that,  in  general,  h  is  not  constant  along  a  tube  but  may  vary  markedly 
with  x/d,  especially  near  the  tube  entrance;  h  tends  to  decrease  from  a 
high  value  near  the  entrance  and  to  approach  a  constant  value  for  large 
values  of  x/d  (14,23).  The  primary  cause  of  this  variation  is  that  the 
temperature  profile  across  the  flow  stream  is  not  fully  developed  at  the 
entrance  of  a  heated  section  (even  though  the  velocity  profile  may  be 
developed  at  the  entrance).  Consequently,  as  the  temperature  profile 
becomes  developed  downstream  from  the  entrance,  the  heat  transfer  coefficient 
approaches  an  asymptotic  value.  Although  the  effect  is  not  as  pronounced  in 
the  case  of  turbulent  flow  as  in  the  case  of  laminar  flow,  it  is  plainly 
evident  from  existing  experimental  results  (14). 

The  variation  of  the  heat  transfer  coefficient  or,  alternately,  of 
the  Nusselt  number  in  the  entrance  region  depends  to  some  extent  on  the 
entrance  Reynolds  number.  However,  it  has  been  demonstrated  (13,  18)  that 
the  variation  can  usually  be  represented  by  a  simple  function  of  x/D  for  a 
wide  range  of  Reynolds  numbers.  Two  different  functions  have  been  used  to 
correlate  average  heat  transfer  coefficients  in  the  case  of  a  fully  developed 
velocity  profile  and  a  flat  temperature  profile  at  the  start  of  heating. 

The  first  (7,10,13)  assumgs^that  the  average  Nusselt  number  (up  to  the 
length  x)  varies  as  until  x/D  =  60;  beyond  this  value  the  Nxisselt 


*This  discussion  supersedes  Section  3.1  of  Technical  Report  247-1  (2). 
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number  is  essentially  constant  with  x/d.  For  values  of  x/D  greater  than  60, 
the  factor  (x/d)"®*^  is  replaced  by  the  constant  (60)"^*^.  The  second  method 
of  correlation  (18,  22)  uses  the  factor 


to  express  the  variation  of  average  Niisselt  number  with  axial  position.  It 
should  be  noted  that,  at  very  large  x/d,  the  two  functions  differ  by  a 
factor  of  (60)”^*^  =  0.664. 

Since  the  local  heat  transfer  coefficient  h,  at  a  point  axray  from  the 
entrance,  is  alijays  smaller  than  the  average  coefficient  h^^  up  to  that 
point,  the  entrance  length  for  the  former  is  considerably  shorter  than  for 
the  latter  and  the  correlations  given  above  can  not  be  rigorously  valid  for 
local  considerations.  However,  no  entrance  length  correlations  are  available 
for  local  heat  transfer  coefficients  and  Nusselt  numbers;  either  local  values 
are  reported  only  at  x/d  greater  than  60  (the  so  called  "fully  developed" 
values)  or  one  of  the  average  correlations  is  used  without  comment. 

3.2  CORRELATION  OF  HEAT  TRANSFER  DATA* 

At  moderate  temperature  differences,  it  has  been  customary  (13)  to 
correlate  heat  transfer  data  for  fluids  flowing  turbulently  inside  smooth 
tubes  by  the  equation 


11 


However,  the  investigation  of  Drexel  and  McAdams  (13)  suggests  that  equation 
11  should  be  modified  slightly  for  air,  and  presumably  other  gases,  to  read 


12 


Equations  used  to  correlate  data  under  conditions  of  large  temperature 
differences  should  be  compatible  with  these  relations .  Compatibility  can 


Supersedes  Section  3.3  of  Technical  Report  247-1  (2). 
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be  achieved,  as  stated  previously,  either  by  evaluating  all  fluid  properties 
at  a  "film”  temperature  between  the  bulk  and  the  surface  temperature  or  by 
introducing  a  temperature  ratio  parameter, 

Humble,  Lowdermilk  and  Desmon  (10),  in  reporting  experimentally 
measured  average  heat  transfer  coefficients  for  air  flowing  in  smooth  heated 
tubes,  give  the  relation 


=  <Vf,av> 


-0.1 


13 


vdiere  (60)"*^*^  was  used  for  values  of  l/D  greater  than  60 .  The  subscriot  f 
denotes  that  flxiid  properties,  including  the  density,  were  evaluated  at  a 
"film"  temperature  equal  to  the  average  of  the  sxirface  and  gas  bulk 
temperatures.*  Equation  13  fits  the  experimental  data  very  well,  especially 
for  Reynolds  numbers  of  10,000  and  above.  The  same  relation  for  average 
heat  transfer  coefficients  was  used  by  Durham,  Neal  and  Newman  (7)  for 
correlating  data  obtained  by  heating  helium  flowing  in  smooth  tubes.  Taylor 
and  Kirchgessner  (l8)  report  a  slightly  different  form  of  the  equation  for 


helium  floijing  in  smooth  tubes: 

’^Nu,f,av  ~  ^\e,f,av^  ^^Pr,f,av^ 


/t  '-0*7 


14 


To  correlate  local  heat  transfer  coefficients  in  the  fully  developed 
region  (x/D  greater  than  60),  Taylor  and  Kirchgessner  used 


^Nu,f  ^Re,f 


0.8 


N 


0.4 


Pr,f 
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*Note  that  the  "modified" Reynolds  number  is  defined  ty 

T 

M  =  b 

\e,f  vj  T^ 


and,  for  a  perfect  gas, 


Dp^u 

^Re,f  ~ 
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which  fits  their  data  within  +  38^.  Lapides  and  Goldstein  (12),  using  a 
corrected  value  for  the  thermal  conductivity  of  air,  replotted  the  fully 
developed  local  data  obtained  by  the  NACA  and  reported  by  Deissler  and  Eian 
(6).  The  best  fit  was  given,  within  approximately  +  105?,  by 


Nu, 


f  =  0*0^5  Nj^e,f 


0.8 


N 


0.4 


Pr,f 
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Durham,  Neal  and  Newman  report,  for  helium  flowing  in  smooth  tubes,  the 
eqiiation 


NNu,f  = 


0.8 


N 


0.4 


Pr,f 


17 


which  reduces  to 


%u,f  =  %e.f 


0.8 


N 


0.4 


Pr,f 
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for  local  Nusselt  numbers  in  the  fully  developed  region. 

Correlations  of  the  type  above,  i.e.,  those  which  employ  the  evaluation 
of  flirid  properties  at  a  "film"  temperature,  tend  to  obscure  the  significance 
of  the  actual  fluid  properties.  In  turbulent  flow  the  velocity  and  temperature 
are  relatively  constant  across  the  turbulent  core,  and  the  gas  properties, 
consequently,  are  also  almost  constant  across  this  turbulent  core.  I&x 
particular  the  mass  flow  velocity  G,  a  quantity  •tdiich  is  independent  of 
distance  in  the  axial  direction  and  which  is  widely  used,  is  not  directly 
employed  in  the  correlations. 

These  anomalies  may  be  avoided  by  using  a  correlation  in  which  the 
fluid  properties  are  evaluated  at  gas  bulk  temperature  and  into  idiich 
is  introduced  as  a  parameter  (l,l4,17).  For  average  Nusselt  numbers, 

Barnes  (1)  reports  the  equation 


N„  =  0,023 
Nu,av 


N. 


0.8 


Re,av 


N 


0.4 


Pr,av 
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where  the  exponent,  a,  varies  from  gas  to  gas.  His  investigation  gives  the 
following  values; 
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Air  a  =  -0.4 

Helium  a  =  -0.185 

Carbon  Dioxide  a  =  -0.27 


At  the  same  time,  Barnes  states  that  the  data  of  Humble,  et  al.  (10), 
suggest  a  =  -0.55  for  air. 

The  results  of  the  present  tests  are  compared  to  both  types  of  correla¬ 
tions  discussed,  i.e.,  those  using  a  reference  "film"  temperature  and  those 
employing  a  temperature  ratio  parameter. 


3.3  local  friction  FACTORS 

Deissler  (3)  reports  analytically  determined  friction  factors  for 
fully  developed,  turbvilent  flow  through  smooth  tubes  with  heating  for  a 
Prandtl  number  of  one.  The  results  indicate  that  friction  factors,  for  the 
case  of  gas  heating,  fall  below  the  generally  accepted  Karman-Nikuradse 
relation  for  isothermal,  turbulent  flow. 
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when  the  fluid  properties  are  evaluated  at  gas  bulk  temperature.  Further¬ 
more,  it  is  indicated  that  this  trend  is  progressively  more  pronounced  at 
T 

higher  values  of  .  The  analysis  predicts  that  equation  20  will  correlate 

'A'U 

fully  developed  flow  friction  factor  data  if  gas  properties  are  evaluated  at 
a  reference  temperature,  Tq  defined  as 
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Deissler  and  Eian  (6)  extended  the  analysis  descidbed  above  to  the  case 
of  a  Prandtl  number  of  0.73,  and  the  resulting  friction  factors  were  compared 
to  experimentally  determined  values  measured  at  a  position  6  inches  from  the 
exit  of  an  87- inch  test  section  for  the  flow  of  air.  The  results  are 
correlated  by  equation  20  for  Reynolds  numbers  above  30,000  if  the  gas 
properties  are  evalxiated  at  the  reference  temperature,  Tq  for  lower 
Reynolds  numbers,  the  experimental  friction  factors  fall  above  those  pre¬ 
dicted  by  the  equation. 
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3.4  AVERAGE  FRICTION  FACTORS 

A  number  of  authors  (10,18,22)  report  experimentally  determined  average 
friction  factors  for  gases  flovdjig  through  smooth  tubes  with  high  surface  to 
bulk  temperature  ratios.  These  results  generally  indicate  that  if  the  fluid 
properties  are  evaluated  at  gas  bulk  temperature,  the  average  friction 
factors  fall  below  those  predicted  by  equation  20  for  Reynolds  numbers  above 
30,000.  These  results  show  that  if  the  fluid  properties  are  evaluated  at 
the  average  "film"  temperature,  the  friction  factors  are  fairly  well  pre¬ 
dicted  ly  equation  20  for  Reynolds  numbers  above  30,000.  At  lovjer  Reynolds 
numbers,  the  friction  factors  tend  to  be  progressively  h^her  than  those 
predicted  by  the  Kanu^-Nikuradsa  relation  as  the  ratio  is  increased, 
idiether  the  fluid  properties  are  evaluated  at  bulk  temperature  or  at  a  film 
temperature . 
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4.  lXPERim}^i:AL  APPARATUS 


Heat  transfer  and  pressure  drop  neasurerients  are  made  for  a  gas  flowing 
through  a  l/8-inch  i.d.  by  0.022-inch  wall,  round,  inconel  tube  x»dth  a 
heated  length  of  20  inches  and  an  unheated  entrance  length  of  2  l/8  inches. 
Thermocouples  and  pressure  taps  are  attached  to  this  test  section  tube  at 
several  axial  positions. 

The  apparatus  described  in  Technical  Report  Mo.  247-1(2)  x-jas  modified 
for  the  present  experimental  program  as  follows: 

1.  A  support  stand  for  the  test  section  added. 

2.  The  electrode  arrangement  \wls  altered. 

3.  Thermocouple  attachment  and  test  section  insulation  vrere  changed. 

4.  The  entrance  gas  mixer  ^;as  placed  inside  the  pressure  vessel  at 
the  inlet  to  the  test  section. 

5.  Pressure  taps  were  added  to  the  test  section,  and  a  pressure  drop 
measuring  system  vra-s  installed. 

4.1  TEST  SECTION  MODIFICATIONS 

Figure  1  shows  the  test  section  assembly  less  thermocouples,  insulation, 
and  pressure  vessel.  The  support  stand  consists  of  three  1/2- inch  inconel 
tubes  screvjed  into  the  lower  flange,  two  reinforcing  brackets,  and  a  plate 
which  supports  the  upper  electrode  assembly.  A  graphite-packed  sliding 
electrode  (2)  is  welded  to  an  inconel  plate  to  which  two  3/8-inch  diameter, 
nickel  bus  bars  are  welded.  This  plate  is  bolted  to,  but  insulated  from, 
the  upper  plate  of  the  support  stand. 

A  3/8-inch  o.d.  by  l/8-inch  i.d.  tube  is  brazed  onto  the  upper  end  of 
the  test  section.  The  tube  is  free  to  move  vertically  in  the  sliding 
electrode  as  thermal  expansion  of  the  test  section  dictates.  The  support 
stand  protects  the  test  section  while  the  pressure  vessel  is  being  lowered 
into  place  and  serves  as  a  frame  for  the  attachment  of  the  pressure  taps  and 
thermocouples , 

Chromel-alximel  thermocouples  are  welded  to  the  test  section  at  nine 
aydal  locations  and  are  supported  by  one-inch  diameter,  alumina  segments 
(Figure  2),  These  segments  are  clamped  to  the  test  section  independently; 
therefore,  they  move  wdth  the  test  section  as  it  experiences  thermal  expan¬ 
sion,  and  the  thermocouples  and  pressure  taps  are  not  torn  from  the  test 
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section.  Quartz  fiber  is  placed  in  the  annulus  between  the  segmented  insula¬ 
tion  and  the  pressure  vessel  to  reduce  the  heat  loss. 

Pressure  taps  are  located  along  the  test  section  at  seven  axial  positions 
so  that  the  static  pressure  variation  can  be  measured  as  the  gas  flows  through 
the  test  section.  These  pressure  taps  are  l/l6-inch  o.d.  inconel  tubes 
brazed  to  the  test  section  with  0.050-inch  holes  drilled  through  the  test 
section  wall.  All  burrs  inside  the  test  section  resulting  from  the  drilling 
process  were  removed. 

A  gas  mixer  is  located  inside  the  pressure  vessel  so  that  the  gas 
temperature  can  be  measured  immediately  upstream  from  the  test  section 
entrance . 

4.2  PRESSURE- DROP  MEASURING  SYSTEM 

A  bank  of  twelve  50-inch  manometers,  a  60-inch  inclined  manometer,  and 
two  0-1200  psi  Bourdon  tube  pressure  gauges  are  used  to  measure  the  static 

pressure  distribution  along  the  test  section. 

The  bank  of  manometers  is  so  arranged  that  the  static  pressure  drop 
between  successive  pressure  taps  and  between  the  last  pressure  tap  and  the 
pressure  vessel  can  be  measured  by  50-inch  water  manometers;  or  the 
successive  pressure  drops  between  the  last  five  taps  and  the  pressure  vessel 
can  be  measured  by  50-inch  mercury  manometers.  An  additional  manifold  and 
valve  arrangement  permit  the  measurement  of  pressure  drop  between  any  two 
taps,  or  any  tap  and  the  pressure  vessel,  by  use  of  either  the  60-inch 
inclined-tube  water  manometer;  a  50-inch  mercury  manometer;  or  two  Bourdon 
tube  gauges.  With  this  arrangement,  pressure  drops  from  0,01  inch  of  water 
to  several  hundred  psi  can  be  measured. 
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5.  METHOD  OF  TESTING 


A  series  of  experiments  has  been  performed  with  nitrogen  to  deteimine 
local  heat  transfer  coefficients  and  local  and  average  friction  factors  for 
gas  flowing  through  tubes  of  circular  cross-section. 

The  method  of  testing  was  the  same  as  reported  in  Technical  Report  No. 
247-1(2)  except  as  described  below.  A  number  of  the  important  expeidmentally 
determined  quantities  are  shown  in  FigiU'es  8  and  9. 

5.1  RANGE  OF  VARIABLES 

Heat  transfer  and  pressure  drop  data  were  obtained  over  the  ranges  of 
important  variables  listed  in  the  following  table: 


RANGE  OF  VARIABLES 


1. 

Flow  rate 

1.24  to  32.6  Ib/hr 

2. 

Inlet  gas  temperature 

approximately  70°  F 

3. 

Maximum  gas  temperature 

784°  F 

4. 

Maximum  surface  temperature 

1455°  F 

5. 

Maximum  T^T-^^ 

2.08 

6. 

Bulk  Re3rnolds  number 

1830  to  90,000 

7. 

Heat  fltix 

733  to  138,000  Btu/(hr)(sq  ft) 

8. 

Axial  location  parameter,  ^ 

2  to  127 

5.2  THERfOCOUPLE  CALIBRATION 

The  test  section  ten^jerature  at  several  locations  was  measured  by 
chromel-alumel  thermocouples  welded  to  the  outside  of  the  test  section. 

Part  of  the  signal  from  a  given  thermocouple  was  caxised  by  the  heating 
ciirrent  flow  through  the  test  section;  i.e.,  the  two  wires  of  the  thermo¬ 
couple  were  not  at  exactly  the  same  potential  on  the  test  section.  This 
effect  was  minimi ted  by  careful  attachment  of  the  thermocouples ^  but 
could  not  be  completely  elimjjiated.  Each  tb.ermocouple  was  calibrated  for 
this  erroneous  sJ^gnal  as  a  function  of  test  section  cufrrent;  a  description 
of  th®  calibration  is  presented  in  Appendix  A,  and  a  sample  resulting  cali¬ 
bration  curve  is  shown  in  Figure  10.  A  typical  test  section  temperatxire 
distribution  is  shown  in  Figure  8.  Inside  wall  temperatures  were  calculated 
as  described  in  Appendix  B, 
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5.3  HEAT  LOSS  CALIBRATION 

The  data  for  heat  loss  calibration  were  taken  the  same  way  as  reported 
In  reference  2;  however,  the  interpretation  of  the  heat  loss  and  subsequent 
application  of  the  results  to  the  experiments  were  refined  as  discussed  in 
Appendix  C. 


^.4  HEAT  yr.UX.  TO  THE  GAS 

The  total  power  input  to  the  test  section  was  measured  as  described  in 
reference  2  except  that  the  voltage  drops  between  the  lower  electrode  and 
the  pressure  taps  were  measured;  these  showed  that  the  voltage  varied  linearly 
along  the  test  section.  The  conv>-rsion  of  electrical  energy  to  thermal 
energy,  consequently,  was  uniform  along  the  test  section.  The  calculation 
of  the  heat  loss  is  described  in  Appendix  C,  and  the  difference  between 
the  heat  generation  per  unit  length  and  the  heat  loss  per  unit  length  is  the 
heat  transferred  to  the  gas  per  unit  length.  A  typical  plot  of  heat  flux  uo 
the  gas  as  a  function  of  |  is  shown  in  Figure  8 . 


5.5  GAS  TEMPERATURE 

The  inlet  gas  temperature  was  measiu*ed  in  a  mixer  inside  the  pressure 
vessd.  just  at  the  test  section  entrance.  Gas  bulk  stagnation  temperature  at 
any  axial  location  was  calculated  by  equating  the  increase  in  enthalpy  of 
the  gas  to  the  integrated  heat  transfer  to  the  gas,  up  to  that  point,  and 
then  using  an  enthalpy  curve  (Figure  11)  to  detenaine  the  temperature. 

The  static  enthalpy,  H,  at  each  axial  location  was  then  calculated  from 

the  energy  equation, 


u 


where  H  is  the  stagnation  enthalpy  (at  the  known  stagnation  temperature) 
and  u  is  the  average  velocity  at  the  location  as  detennined  from  the  con¬ 
tinuity  equation,  w  =  pAu,and  the  equation  of  state  for  nitrogen,  P  -  P(P,T). 
The  static  temperature  was  read  from  a  plot  of  H  vs  T  with  pressure  as  a 
paimeter.  Since  the  density  is  a  function  of  static  temperature,  the 
density  and  average  velocity  could  not  be  found  until  the  static  temperature 
IAS  known,  -nierefore,  it  was  necessary  to  use  an  iterative  procedure;  the 
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static  temperature  vras  estimated,  density  and  velocity  determined,  and,  from 
the  static  enthalpy,  a  corrected  temperature  was  found,,  The  calculation 
converged  very  rapidly,  only  two  or  three  cycles  being  required  at  each 
location . 

In  addition  the  adiabatic  wall  temperature  was  calculated  hy  the  equation, 
Ta^  =  T  +JR(T^  -  T)  =  T  +  -  T) 

given  in  reference  15  where  the  exponent,  a,  is  l/3  for  turbulent  flow  and 
1/2  for  laminar  flow.  The  turbulent  value  was  used  for  all  runs;  however, 
no  significant  error  was  hereby  introduced  in  the  low  Reynolds  number  runs 
as  the  maximum  difference  between  static  and  stagnation  temperatures  was 
only  37°  F. 

5.6  PRESSURE  DROP 

The  pressure  drops  between  successive  pressure  taps  and  between  the 
last  tap  and  the  pressure  vessel  \jere  measured  during  each  run.  The  coiri- 
bination  of  manometer  ban!:,  inclined  manometer,  and  Bourdon  tube  gauges 
was  used  in  such  a  manner  as  to  measure  each  pressure  drop  with  maximum 
available  sensitivity;  i.e.,  the  inclined  manometer  was  used  to  measure 
pressure  drops  from  0  to  10  inches  of  water,  the  water  manometers  to  measure 
pressure  drops  from  10  to  50  inches  of  water,  the  mercury  manometers  to 
measure  pressure  drops  from  50  inches  of  water  to  50  inches  of  mercury, 
and  the  Bourdon  tube  gauges  to  measure  pressure  drops  greater  than  50 
inches  of  mercury.  Figure  8  shows  a  typical  curve  of  static  pressure 
variation  along  the  test  section  with  the  vessel  pressure  as  the  datum. 

5.7  HEAT  TRANSFER  COEFFICIENTS 

Using  the  calculated  heat  flux  to  the  gas,  q^,,  and  the  temperature 

difference  T  -  T  .  the  local  heat  transfer  coefficient,  h,  was  calcula- 
’  w  aw’ 

ted  from  equation  1, 


A  typical  plot  of  the  heat  transfer  coefficient  for  a  test  is  shown  in 
Figure  9. 
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5.3  FRICTION  FACTORS 

Average  friction  factors  were  calculated  from  equaiio;:is  8  and  9  with 


the  two  stations  being  pressure  taps  1  and  7.  Local  friction  factors  were 

H  r  pu'^  1  -  . 


calculated  from  equations  5  and  6.  The  tern,  ^ 
at  selected  thermocouple  positions  by  plotting 


P  + 


was  evaluated 


P  + 


pu 


as  a 


function  of  x  and  then  graphically  differentiating  the  resulting  curve  at 

the  position  of  interest.  A  typical  plot  of  local  friction  factor  as  a 

function  of  ^  is  shovm  in  Figure  9. 

D 


5,9  GAS  PROPERTIES 

Four  properties  of  nitrogen  were  used  in  the  calculations:  enthalpy, 
viscosity,  thermal  conductivity,  and  Prandtl  number.  These  properties  were 
evaluated  at  the  associated  static  bulk  or  film  temperature  except  in 
determining  the  stagnation  enthalpy  xdiere,  by  definition,  the  stagnation 
temperature  must  be  used.  Figures  11  through  Ik  in  Appendix  E  show  these 
properties*  variation  >7ith  temperature,  and  the  important  references  are 
listed. 
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6.  EXPERIi'SNTAL  RESULTS 


The  results  of  the  present  experimental  investigation  are  graphically 
summarized  in  Figures  3,  and  5  and  are  discussed  in  detail  below.  The 
heat  transfer  results  are  presented  in  two  ways.  In  the  first  (Figure  3), 
the  gas  properties  are  evaluated  at  the  "film”  temperature,  as  defined  in 
Section  3.  In  the  second  method  (Figure  4),  the  bulk  static  temperature  of 
the  gas  is  used  for  property  evaluation  and  the  temperature  ratio,  T  Jl,, 
is  introduced  as  a  parameter.  Fully  developed  local  friction  factors  are 
presented  in  Figures  5  and  6  with  gas  properties  evaluated  at  bulk  and  film 
temperature  respectively.  Figiire  7  is  a  plot  of  average  friction  factors. 

In  addition.  Figures  8  and  9  show  representative  axial  distributions  of  the 
more  important  experimental  variables. 

Because  of  their  large  uncertainties,  the  data  obtained  with  an  earlier 
test  section  and  presented  in  Technical  Report  247-1(2)  are  not  included  in 
the  present  report.  In  general,  the  heat  transfer  correlations  used  in 
report  247-1  lie  somewhat  above  those  discussed  below  in  Sections  6.1  and 
6.2. 


6.1  HEAT  TRANSFER  CORRELATION  BASED  ON  FILM  TEMPERATURE 

The  local  heat  transfer  data  are  correlated,  as  shoxm  in  Figure  3,  tiy 
plotting  the  quantity 


N 


Nuj;_ 


N 


0.4 


Pr,f 


1  + 


as  a  function  of  the  film  Reynolds  nxmber.  The  best  fit  to  the  data  was 
given  by 


%u,f  ^Re,f 


0.8 


N 


0.4 


Pr,f 
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which  is  shown  as  a  solid  line.  The  two  dashed  lines  indicate  10  per  cent 
deviation  from  this  equation.  The  bulk  of  the  data  falls  within  the  +  10 
per  cent  interval,  the  exceptions  being  data  points  from  the  two  thermocouples 
at  ^  smaller  than  6,  where  axial  conduction  heat  losses  are  large. 


Modified  Reynolds  Number,  NRe,t 


Figure  3 


Heat  Transfer  Correlation  Based  on  Local  Values, 

Fluid  Properties  Evaluated  at  Reference  Temperature,  T^ 


Only  those  iruns  with  film  Reynolds  numbers  greater  than  8000  are  plotted. 
However,  the  tendency  of  the  data  to  fall  below  the  turbulent  flow  correlation 
at  low  Reynolds  numbers  is  already  evident.  This  trend  is  to  be  expected 
as  the  flow  may  be  transitory  between  the  laminar  and  turbulent  regimes  of 
flow. 

The  experimental  curve  (equation  22)  is  slightly  lower  than  those 
reported  by  other  investigators,  but  it  agrees,  within  the  spread  of  the  data, 
with  eqxuition  l6  for  air  (12)  and  equation  15  for  helium  (l8).  Equations  17 
and  18  taken  from  Durham,  Neal  and  Newman  (?)  lie  considerably  above  the 
present  experimental  correlation  and  others  cited. 


6.2  HEAT  TRANSFER  CORRELATION  BASED  ON  T^T^^  PARAMETER 

The  local  heat  transfer  data  are  correlated  in  terms  of  the  parameters 
^Nu’  %e’  ^Pr’  i  V'^b  ^  Plotting  the  quantity 


“•“E  Hi) -'’•!] 


as  a  function  of  the  Reynolds  number  in  Figure  4.  All  data  obtained  in  the 
present  investigation  are  included  in  this  figure,  including  those  near  the 
entrance  where  x/d  is  2.  The  solid  line  represents 


the  best  fit  to  the  data  with  Reynolds  numbers  greater  than  10,000,  and  the 
dashed  lines  indicate  10  per  cent  deviation  from  this  equation.  Again,  the 
data  with  Reynolds  number  greater  than  10,000  are  correlated  within  +  10  per 
cent,  while  the  two  data  points  on  each  run  where  ^  is  less  than  6  tend 
toward  greater  scatter. 

The  range  of  values  of  T^T^  in  the  present  experimental  results 
(maximum  of  2.1)  was  insiifficient  to  check  the  generality  of  the  exponent 
-0.5  appearing  on  T in  equation  23.  However,  the  data  do  indicate  that 
the  value  is  somewhat  greater  than  the  -0.4  reported  by  Barnes  (l)  for  air. 
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Reynolds  Number,  Nrq 


Figure  4 


Heat  Transfer  Correlation  Based  on  Local  Values 
Fluid  Properties  Evaluated  at  Bulk  Temperature 


6.3  THERMAL  ENTRANCE  REGION 

Both  of  the  entrance  parameters  discussed  in  Section  3*1, 


and  1  + 


vrere  investigated  as  means  of  correlating  the  experimental  data.  There 

yi. 

is  no  noticeable  difference  between  the  two  parameters  for  j:  greater  than  6, 

X  /x\-0  1  ^ 

At  smaller  the  two  diverge  with  (-j  *  falling  well  helo\i 


On  the  other  hand,  the  data  points  from  the  first  two  thermo¬ 


couples  =  1.7^  and  5.7^)  lie  slight]^  above  the  values  predicted  using 


the  latter  parameter.  Consequently, 
heat  transfer  correlation  equations. 


was  chosen  for  the  two 


6  A  FRICTION  FACTORS 

The  experimentally  determined  friction  factors  are  plotted  vs  Reynolds 

number  on  Figures  5,  6,  and  7  with  Karman-Nikuradse  relation  for  isothermal, 

turbulent  flow  through  smooth  tubes  and  the  laminar  flow  relation, 

f  =  ,  included  for  reference.  The  results  of  a  series  of  tests  per- 

Re 

formed  to  measure  average  friction  factors  with  no  gas  heating  are  shown 
in  Figure  7  and  indicate  that  the  test  section  may  be  considered  as  a  smooth 
tube  for  friction  factor  measurements. 

Local  friction  factors  as  defined  by  equation  5  were  determined  at 
each  thermocouple  location,  and  Figure  9  shows  a  plot  of  friction  factor 
variation  along  the  test  section  for  a  typical  run.  As  can  be  seeh  in  the 
figure,  the  friction  factor  increases  in  the  direction  of  flow  in  the 
thermal  entrance  region  and  then  decreases  in  the  direction  of  flow  toward 
the  exit.  The  large  variation  in  the  friction  factor  along  the  test  section 
clearly  illustrates  that  the  pressure  drop  is  significantly  affected  hy  the 
heat  transfer  process,  in  addition  to  the  effect  of  a  simple  change  in  the 
momentum  of  the  gas.  It  is  also  evident  that  the  flow  cannot  be  considered 
"fxilly  developed"  in  a  rigorous  sense,  even  at  large  values  of  ^  ,  as  far 
as  the  friction  factors  are  concerned. 

The  friction  factors  calculated  at  thermocouple  number  7  (^  =  9^), 
which  is  as  near  the  exit  as  end  effects  can  be  neglected  for  all  runs,  are 
plotted  vs  Reynolds  number  with  gas  properties  calculated  at  bulk  temperature 
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Figure  5  Friction  Factor  Correlation  Based  on  Local  Values, 
Fluid  Properties  Evaluated  at  Bulk  Temperature 


Present  Tests 
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Fluid  Properties  Evaluated  at  Bulk  Temperature 


Figure  8  Heat  Flux,  Temperature,  and  Pressure 

Distributions  for  a  Typical  Test 
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Static  Pressure  Above  Datum 


Figure  9 


Distance  Parometer,— 

D 

Heat  Transfer  Coefficient,  Friction  Factor 
Nusselt  Number,  Reynolds  Number,  and 
Temperature  Ratio  Variation  for  a 
Typical  Test 


in  Figure  5.  For  P^ynolds  numbers  above  10,000,  the  results  are  substantially 
in  agreement  with  those  predicted  Deissler  (5)  For  fully  developed  flow 
in  that  the  friction  factors  are  lovjer  than  would  be  calculated  by  equation 
20.  At  the  lower  Reynolds  numbers,  the  friction  factors  are  considerably 
higher  than  those  predicted  by  equation  20  or  the  laminar  flow  relation. 

Figure  6  shows  the  same  points  but  id.th  the  fluid  properties  evaluated  at 
a  reference  temperature,  T^;  the  results  are  correlated  by  equation  20  for 
Reynolds  numbers  above  10,000  and  are  above  the  correlation  for  lower  Reynolds 
numbers.  These  same  trends  are  observed  in  the  results  of  Deissler  and 
Eian  (6). 

Figure  7  shows  the  average  friction  factor  as  a  function  of  the  Reynolds 
number  with  all  fluid  properties  evalixated  at  the  average  gas  bulk  temperature. 
The  results  of  Taylor  and  Kirchgessner  (18)  are  also  plotted  on  the  figure. 

For  Reynolds  numbers  above  20,000,  the  friction  factors  are  slightly  below 
the  Karman-Nikuradse  relation;  and  for  lower  Reynolds  numbers,  the  friction 
factors  are  considerably  above  this  relation  as  well  as  above  the  laminar 
relation.  No  attempt  was  made  to  correlate  average  friction  factor  data 
Trri.th  fluid  properties  evaluated  at  an  average  reference  temperature  becaxise 
there  is  no  significant  average  reference  temperat’ore  wiien  the  conditions  are 
changing  so  much  along  the  test  section  (Figures  8  and  9). 

6.5  EXPERIMENTAL  UNCERTAINTf 

The  uncertainties  in  the  reported  experimental  results  were  determined 
from  the  estimated  uncertainty  inteinrals  associated  with  the  various  primary 
measurements.  Estimates  of  the  various  uncertainty  intervals  were  obtained 
from  instrument  specifications  and  calibrations,  observed  measurement 
fluctuations,  and  operator  judgement.  A  detailed  uncertainty  analysis,  after 
the  method  of  reference  11,  is  presented  in  Appendix  D.  The  numerical  values 
of  the  uncertainty  intervals  obtained  are  reasonable  estimates  of  typical 
intervals  for  the  present  tests. 

An  uncertainty  interval  of  +  9.5  P®^  cent  was  found  for  the  heat 
transfer  coefficients  and  is  reliable  except  for  the  results  obtained  near 
the  lower  electrode,  vdiere  the  uncertainty  interval  is  considerably  larger 
because  of  the  significant  axial  conduction  heat  losses. 
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The  uncertainty  interval  in  the  friction  factors  is  +  12  per  cent,  and 
the  controlling  item  is  the  numerical  or  graphical  differentiation  necessaiy 
to  evaluate  the  apparent  wall  shear  stress  from  equation  7 . 
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APPENDIX  A.  THERMOCOUPLE  CALIBRATION 


Thermocouples  in  electrical  contact  idth  c\irrent- carrying  conductors 
will  produce  a  voltage  signal  consisting  of  two  parts:  the  desired  emf  ^ich 
is  approximately  proportional  to  the  temperature  difference  between  the  hot 
and  cold  junctions,  and  an  extraneous  component  ^ich  occurs  because  the  two 
wires  constituting  the  thermocouple  junction  are  not  at  the  same  electrical 
potential  on  the  conductor*  In  order  to  determine  the  temperature  in  such 
a  situation,  the  extraneous  signal  must  be  measured  or  eliminated.  For 
example,  if  the  cxirrent  in  the  conductor  is  alternating,  an  it-C  filter 
circuit  may  be  devised  to  attenuate  this  component. 

In  welding  the  thermocouples  to  the  test  section,  care  was  taken  to  keep 
the  two  wires  at  the  same  electrical  potential  hy  trimming  the  junction  and 
carefully  positioning  it  on  the  test  section  before  welding.  Since  this 
pix)cedure  did  not  completely  eliminate  the  extraneous  signal,  a  method  of 
calibration  was  devised  to  measure  it  accurately  as  a  function  of  test 
section  current  and  temperature. 

The  extraneous  signal  was  measured  by  passing  a  steady  heating  current 
of  up  to  300  amperes  through  the  test  section  while  continuously  recording 
the  total  thermocouple  emf  and  the  heating  current  on  an  oscillograph.  The 
current  flow  was  then  interrupted  suddenly  by  opening  a  300-ampere  circuit 
breaker,  after  which  the  thermocouple  emf  no  longer  included  the  extraneous 
component.  Figure  10a  shows  a  replica  of  a  typical  oscillograph  trace  from 
such  a  test,  from  vdiich  it  may  be  seen  that  the  cooling  of  the  test  section 
after  the  interruption  of  heating  current  is  slow  compared  to  the  rate  of 
change  of  heating  current.  Consequently,  the  sudden  change  in  thermocouple 
emf,  ^ich  may  be  positive  or  negative,  is  entirely  attributable  to  the 
extraneous  pickup. 

The  calibration  procedure  descii.bed  above  was  repeated  several  times 
for  each  thermocouple  with  different  values  of  heating  current  and  wall 
temperature.  Figure  10b  shows  the  extraneous  signal  plotted  as  a  function 
of  the  current  for  one  thermocouple  with  wall  temperatures  from  85  to 
1270^F.  From  tests  such  as  these  on  all  the  thermocouples,  it  is  concluded 
that  the  pickup  from  the  test  section  varies  linearly  with  heating  current 
and  is  independent  of  temperature. 
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Section  Current^  Amperes 


Thermocouple  Signal,  Millivolts 


APPENDIX  B.  INSIDE  SURFACE  TEMPERATURE 


Thermocouples  welded  to  the  outside  of  the  test  section  tube  are  the 
only  means  of  test  section  temperature  measurement.  To  obtain  inside  sur¬ 
face  temperatures,  it  is,  of  course,  necessary  to  correct  these  temperature 
measurements  for  temperature  variation  across  the  wall  of  the  test  section. 
The  method  of  correction  described  below  was  also  used  in  Technical  Report 
24’7-l(2)  with  the  exception  that,  in  Report  24’7-l,  the  axial  conduction 
heat  loss  was  neglected. 

Since  the  test  section  assembly  is  symmetric  with  respect  to  the  angular 
coordinate  at  a  given  axial  location,  the  conduction  equation  for  the  test 
section  is  independent  of  the  angular  coordinate.  If  axial  conduction  is 
neglected  (this  term  is  small  except  near  the  upper  and  lower  electrodes), 
the  conduction  eqiiation  will  be  independent  of  the  axial  coordinate  and  the 
equation  may  be  written  as 


1  ^ 
r  dr 


B-1 


where  T  is  the  local  temperature  in  the  tube  wall,  r  is  the  radial  coordinate, 
is  the  thermal  conductivity  of  the  tube  material,  and  q'**  is  the  heat 
generated  per  unit  volume  of  tube  material,  by  the  electrical  energy  dissi¬ 
pated  in  the  test  section. 

If  uniform  heat  generation  is  assumed  and  the  following  boundary  condi¬ 
tions  are  used. 


(a)  q^  =  q^  (from  heat  loss  calibration) 

(b)  T^  =  T  (from  thermocouples) 


the  solution  to  equation  B-1  is 


T  =  T  _  + 
w  wo 
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where  r  and  r  are  the  inside  and  outside  radii  of  the  test  section, 
w  wo 

respectively. 
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The  effect  of  axial  conduction  near  the  electrodes  vras  approximated  by- 
lumping  the  axial  heat  loss  with  the  radial  loss.  That  is,  the  ff.rst  boundary 
condition  was  modified  as  follows: 

(a*)  =  q;  + 


Correspondingly,  equation  ^2  was  changed  to 


1 1 1 


T  = 
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wo 


WO 


-  rf)  - 


I  t  I 
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■«o'  ^  "  sf  \ 


B-2a 


This  equation  was  used  to  calculate  the  inside  surface  temperature,  T^, 
frOTi  the  expeidmental  da-ta. 
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APPENDIX  C.  HEAT  LOSS  CALIBRATION 

At  steady  state,  the  energy  added  to  the  gas  per  vinit  length  of  test 
section,  q^,  must  equal  the  total  energy  generated  in  the  test  section  per 
unit  length,  q^,  minus  the  two  heat  losses:  radial  loss  through  the  s\irround- 

ing  insulation,  q®.,  and  axial  loss  ty  conduction  along  the  inconel  test 
section,  q* .  at  any  axial  position  along  the  test  section, 

q®.  "  q.p  “  (ql  +  qp 

W  i  Ad 

The  iiidependent  determinatior!  of  the  t^-jo  heat  losses  is  the  subject  of  this 
appen&lx, 

Xt  Kas  found  convenient  to  consider  the  test  section  tube  in  two  parts, 
a  section  close  to  the  lov^er  electrode  (x/D  <  10),  \diere  axial  conduction  is 
of  great  importance,  and  s  central  section  (x/D  >  10)  w^here  axi.a].  loss  is 
much  less  than  radial  loss..  The  part  of  the  tube  near  the  upper  electrode 
ms  not  considered  as  no  data  were  obtained  closer  than  fomr  inches  to  this 
end  of  the  test  section, 

C.l  CENTRAL  REGION.  x/P  GREATER  THAN  10 

Calibration  runs  were  made  in  which  power  was  supplied  to  the  test 
section  with  stagnant  gas  inside  (i.e.,  no  flow).  For  these  runs,  all  enei^ 
generated  in  the  test  section  must  be  dissipated  either  radially  through  the 
insulation  or  axially  through  the  electrodes. 


ql  =  q*  +  q* 
^T  ^r  ^a 


The  corrected  wall  temperature  profiles  were  carefully  plotted  and  the 

second  axial  derivative,  _  ,  at  each  thermocouple  location,  was  graphically 

dx 

calculated.  Since  there  was  only  a  small  variation  in  temperature  across 
the  tube  wall  at  any  location,  the  axial  heat  loss  was  estimated  as 


q!  =  -k.  A, 
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where  is  the  thermal  conductivity  of  the  inconel  tubs  and  is  the  cross 
sectional  area  of  the  tube.  In  all  cases,  oj  was  qui.te  small  compared  to 

d 

ql,  of  the  order  of  10  per  cent  or  less. 

The  radial  heat  loss  was  then  calculated  from  equation  C~2,  and  a 
secondary  quantity,  U(T^^,x),  defined  by 

o'  =  U(T  x)  r T  -  T  ,1  C-4 

was  plotted  as  a  function  of  >all  temperature  at  each  thermocouple  location. 

Thus,  for  a  test  run, the  radial  heat  loss  at  each  thermocouple  position 

could  be  calculated  from  the  measured  wall  temperature  and  the  calibration 

curve  of  U(T  x) ,  The  axial  heat  loss  could  be  determined  in  the  same 
wo’ 

manner  as  for  the  calibration  runs. 


C . 2  ENTRANCE  REGION.  x/P  LESS  THAN  10 

Axial  heat  loss  is  of  primary  importance  in  the  entrance  region.  At 
thermocouple  no.  2  (x/D  =  5.7^),  the  axial  loss  is  approximately  equal  to 
the  radial  loss,  and  at  thermocouple  no.  1  (x/D  =  1.7^)  the  axial  loss  is 
more  than  twice  the  radial  loss . 

Since,  in  this  region,  the  accuracy  with  which  the  second  derivative 
could  be  determined  was  poor,  an  equation  of  the  form: 

T  (x)  =  B[1  -  e"-^]  C-5 

WO  ^ 


vras  fitted  to  each  temperature  profile  in  the  entrance  region  ♦ 
axial  derivative  was  then  calculated  as 


ax^ 


e 


The  second 


C-6 


and  the  axial  heat  loss  as 


q'  =  C-7 

B.  C»  U 

A  semi-- analytical  method  was  used  to  determine  the  radial  heat  loss,, 
The  insulation  surrounding  the  test  section  was  considered  as  a  cylinder  of 


~  36  - 


length  L  with  an  inside  radius  equal  to  that  of  the  test  section,  r  ,  and 
an  outside  radius  extending  to  infinity.  Assuming  constant  thermal  con¬ 
ductivity,  the  differential  equation  for  the  steady  state  temperature 
dstribution  in  the  insulation  x>jas 


5^T  ,  1  5T 

ar^ 


C-8 


The  following  boundary  conditions,  approximating  the  physical  situation  near 
the  lower  electrode  and  test  section,  were  set: 

1.  At  the  lower  electrode  (x  =  0),  the  temperature  is  constant  and 
equal  to  the  experimentally  measured  lower  electrode  temperature. 

2.  At  large  distances  from  the  test  section  (r -*•»),  the  insulation 
temperature  approaches  that  of  the  lo^rer  electrode.  Actually,  it  approaches 
an  ambient  value;  however,  this  assumption  greatly  simplifies  the  solution 
and  introduces  little  error  at  the  test  section  near  the  lower  electrode 
idiere  the  solution  is  applied. 

3.  ^t  the  inside  surface  of  the  qylinder  (r  =  r^),  the  temperature 
is  that  of  the  test  section  as  given  by  the  expression 

T  =  B[1  -  e"^] 
wo  '•  •' 


4.  At  X  L,  0* 

VH.th  these  boundary  conditions,  the  solution  of  equation  C-8  is  a 
Foinrier  sine  series  of  the  form 


DO 

T(r,x)  =  ^  \  (fe  ^o  (2L 


C-9 


n=o 


where  ^  is  the  zero  order,  modified  Bessel  function  of  the  second 

kind.  Since  the  radial  heat  loss  at  the  test  section  must  equal  the  radial 
heat  flow  into  the  insulation, 


T 

‘^r  ■  “^'eff  ^^wo^  ^ 


r=r 


wo 


n=o 
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idiere  Kj^  is  the  first  order,  modified  Bessel  function  of  the  second 

kind. 

Hie  effective  thenioLl  conductivity  of  the  insulstion,  determined 

from  the  calibration  (no  flow)  runs  by  solving  equation  C-10  using  the  first 
seven  terms  of  the  series.  (Actually,  the  first  four  terms  were  used  as  all 
terms  with  n  even  are  zero.)  Then,  with  known,  equation  C-10  was  used 

to  calculate  rada.al  heat  losses  at  the  first  two  thermocouples  for  the 
experimental  runs. 
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APPENDIX  D.  UNCERTAINTr  ANALYSIS 


The  method  and  terminology  in  the  following  analysis  follow  those  of 
ln^wA  aad  McCLlntock  (U);  the  rationale  of  the  procedure  is  described  in 
Appendix  I  of  Technical  Report  2^7-1  (2). 


It  is  necessary  to  examine  separately  the  variables  in  the  equation 


defining  the  heat  transfer  coefficient, 


D-1 


Let  6  Q  represent  the  uncertainly  interval  of  a  variable  G;  then  by  the 
the  method  of  ref ( 

h 

and  substituting  from  equation  D-1, 


where  q^  is  the  total  power  supplied  per  unit  length  of  test  section,  q^  is 
the  radial  heat  loss  per  unit  time  per  unit  length,  q^  is  the  axial  heat 
loss  per  unit  time  per  unit  length,  and  A^  is  the  inside  svirface  area  per 
unit  length.  The  uncertainty  in  the  area  A^  can  be  neglected  in  comparison 
with  the  uncertainties  in  the  energy  terms.  Thus,  the  percentage  uncertainty 
in  the  heat  flux  to  the  gas  can  be  expressed  as 
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Since  the  values  of'  q*  were  small  compared  to  q'  except  near  the  ends  of  the 
heated  section,  the  last  term  can  be  neglected  over  most  of  the  test  section. 
However,  near  the  lower  electrode  (small  x/d),  the  magnitude  of  q^,  coupled 
with  its  large  uncertainty,  resulted  in  a  very  high  heat  flux  uncertainty, 

K 

q"  * 

^  The  uncertainty  in  the  stagnation  temperature  of  the  gas  can,  be  calcu¬ 
lated  as  described  in  Technical  Report  No.  347-1  (2).  The  final  expression 
for  the  absolute  uncertainty  in  the  stagnation  temperature  is 


where  T,  is  the  temperature  of  the  gas  entering  the  test  section  and  q^  is 
the  average  heat  transfer  to  the  gas  per  unit  length.  The  absolute  mcertain- 
ties  of  the  static  and  adiabatic  wall  temperature  can  be  assumed  equal  to 
that  of  the  stagnation  temperature.  The  absolute  uncertainties,  introduced 
apart  from  that  in  the  stagnation  temperature,  are  negligible,  since  the 
three  temperatures,  T,^^  T^  and  differ  by  less  than  40°F. 

As  in  reference  2’  the  inside  surface  temperature  of  the  test  section 
may  be  calculated  by  the  equation 

T  =  T  +  correction  term 
w  wo 

The  correction  terra  is  small,  so  a  large  uncertainty  in  it  would  not 

substantially  affect  the  uncertainty  in  T^;  therefore,  the  xmcertainty  in 

T  xdLll  be  taken  as  the  same  as  that  in  T  .  o 

w  wo 

FRICTION  FACTORS 

The  defining  equation  for  the  friction  factor  can  be  written  as 


2gc  ^^c 


Thus,  the  percentage  uncertainty  in  the  friction  factor  may  be  written  as 
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The  unceirtainty  in  "thn  mass  Plow  velocity  G  is  due  to  the  uncei*ta.5jnly  in  the 
flow  rate  v;  hence  ^  .  Using  the  perfect  gas  relationship 

’  u  V 

P 

^  RT 

the  uncertainty  in  the  density  can  be  related  to  the  uncertainties  in  the 
measured  pressure  and  temperature: 


Special  consideration  must  be  given  to  the  trai-l  shear  stress  uncertainty , 
particularly  for  local  friction  factors .  The  wall  shear  stress  is  given  by 

T 

W 

The  uncertainty  in  the  wall  shear  stress  includes  the  uncertainties  in  P,  G, 
p,  and  the  operation  of  graphically  determining  the  derivative.  A  reasonable 
estimate  of  the  uncertainty  interval, 5  which  included  these  uncertainties, 
was  made. 

MUMERICAL  EVALUATION  OF  UIv^GERTAINTY  INTERVALS 

The  table  on  the  following  page  contains  a  summary  of  the  estimated 
uncertainties  and  the  subsequent  values  of  uncertainty  interval  in  the  calcu¬ 
lated  results,  h  and  f;  the  necessary  calculations  were  made  as  outlined 
above.  The  uncertainty  intervals  in  the  table  are  to  be  regarded  as  typical 
rather  than  specific  values.  It  should  be  noted  that  the  uncertainty 
presented  for  the  local  heat  transfer  coefficient  is  not  valid  in  a  region 
close  to  the  lower  electrode  <  10^  .  Tnere,  the  large  art.sl.  conduction 
heat  loss,  vjith  its  high  uncertainty,  causes  a  greatly  increased  imcertainty 
in  the  heat  transfer  coefficient. 
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Sumary  of  Uncertainty  Analysis 


Defined  Uncertain'ty  Estimated  IMcertainty 

Variable  _ Interval  Interval _ 


±  25S 

±  5f' 

^a 

^4 

±  50^ 

WjG 

bjr  _  6jG 

w  Cj 

±  1^’ 

6  T 

±  5^'f  ;  ± 

e 

V®  ‘  \e 

±  3'^P  ;  ±  ^o6i 

^b,o 

”b,o 

\o 

+  4.3^^ 

T 

aw 

”aw  »Vo  "aw 

T  T,  ’  T  -  T 

aw  b,o  w  aw 

+  4.35s  ;  ±  QM  * 

p 

bp 

P 

+  25S 

i^p 

B(^P)  ; 

±  .025p^i.  !  i  0.5^ 

Calculated  Values 
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+  3M 

cq; 

±  3.^ 

h 

h 

+  9.5^ 

P 

P 

±k.8i 

St 

T 

V 

w 

"w 

±  11^ 

f 

f 

±  1256  ' 

* 

Calculated  Values 
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APPENDIX  E.  PROPERTIES  OF  NITROGEN 

Tlie  major  reference  for  the  properties  used  in  this  report  is  NBS 
circxilar  564  by  Hilsenrath,  et  al.  (9),  a  critics!  sumriiary  of  gas  property 
data  reported  as  of  the  year  195^.  This  source  has  been  supplemented  by 
more  recent  work,  primarily  the  thermal  conductivity  measurements  of  Vines 

(a). 

1.  Enthalpy 

Hilsenrath,  et  al.,  calculated  enthalpy  over  the  range  of  temperatures 
and  pressures  needed  for  this  report.  Their  values  appear  to  be  as  accurate 
as  any  now  available  and  are  in  excellent  agreement  with  1±ose  recently 
calculated  at  the  Franklin  Institute  (19). 

2 .  Viscosity 

The  calculations  of  Hilsenrath,  et  al.,  are  based  on  the  Leonard— Jones 
6-1,2  model  of  kinetic  theoi^^  and  fit  the  avail.abic  exper.lmental  data,  reported 
at  temperatures  as  liigh  as  2300°F,  within  4;  2  per  cent. 

3 .  Thermal  Conductivity 

As  of  1954,  reliable  experimental  data  for  the.’  thermal  conductivity  of 
nitrogen  extended  only  to  890°F.  Hilsenrath,  et  ai. .,  fit  a  semi empirical 
equation  to  this  data  and  extrapolate  to  1700^F<.  However,  the  extrapolation, 
in  the  range  of  900-l650®F,  relies  heavily  on  the  data  of  Stops  (I6),  which 
are  probably  3  l/2  per  cent  too  low  at  l650*^F.  The  more  recent  investigations 
of  Vines  (21)  and  Timrot  and  Vargaftig  (20)  ejrt.end  to  l650°F  and  are  in  good 
agreement.  Their  data,  and  the  loxr  temperature  data,  are  well  correlated 
by  an  equation  proposed  by  F.  G.  Keyes  and  reported  bjr  Vines.  This  equation 
is  plotted  in  Figure  13. 

4.  Prandtl  Number 

The  Prandtl  number  has  been  calculated  as  uC.„/k  using  the  enclosed 
viscosj.ty  and  thermal  conductix’-ity  curves  and  values  of  specific  heat  talcen 
from  1\!BS  circular  5<54. 
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APPENDIX  F.  4000‘  ’F  EXPERIMENTAL  APPARATUS 


A  new  test  section  assembly  and  pressure  vessel  have  been  designed  and 
constructed  to  measure  local  arid  average  convective  heat  transfer  coefficients 
for  helium  flovdng  through  roiT/d  tubes  mth  gas  temperatures  as  high,  as 
4000 °F.  The  new  components  ax^e  compatible  with  the  present  test  appai-atus 
and  can  be  used  interchangeably  with  the  2500°F  test  section  assembly  and 
pressure  vessel.  Figure  15  shows  the  new  pressure  vessel  prior  to  welding, 
and  Figure  l6  is  a  schematic  drawing  of  the  test  section. 

The  pressxrre  vessel  is  formed  from  a  33-inch  length  of  5  inch  diameter, 
scheduled  80  steel  pipe,  cut  parallel  to  its  axis  to  provide  a  four  inch 
>ri.de  side  opening,  ibe  pipe  is  welded  onto  a  rectangular  flange  and  is  closed 
on  the  ends  by  two  arc-shaped  pieces.  A  one-inch  thick,  flat,  steel  plate 
bolts  to  the  flange,  the  pressure  seal  being  made  with  a  copper  gasket.  A 
row  of  six  view  ports  are  located  on  the  cover  plate  so  that  a  radiaticn 
pyrometer  can  sight  on  the  test  section  tube  inside  the  pressure  vessel. 

The  viewing  windows  are  3/4-ij'Ach  diameter  quartz  glass  cylinders,  each 
3 /4-inch  long. 

Pressure  tap  and  thermocouple  packing  gland  fittings  are  provided  on 
the  vessel's  side.  Electrode  packing  glands  are  located  on  the  txvo  ends  of 
the  vessel,  as  are  fittings  for  the  entrance  and  exit  of  the  test  gas. 

Helium  enters  through  the  lower  end  of  the  pressure  vessel  into  a 
mixing  chamber  viiere  the  biilk  average  ter^Derature  of  the  gas  is  measured 
with  a  chromel-alumel  thermocouple  probe.  The  gas  then  enters  the  test 
section,  a  l/8-inch  diameter  thin-walled  tantalum  tube  with  a  heated  lesigth 
of  20  inches.  An  unheated  entrance  ].ength  of  fom:  inches  is  provided  to 
insure  that  the  gas  velocity  profile  is  fully  developed  by  the  initial  point 
of  heating.  After  it  has  flowed  through  the  test  section,  the  helium  is 
allowed  to  escape  into  the  enclosing  pressure  vessel,  th\js  equalizing  the 
pressures  on  the  inside  and  outside  of  the  test  section  tube. 

The  test  section  txibe  is  resistively  heated  by  passing  an  electi.'ic 
current  through  it.  Current  is  carried  to  the  test  section  through  a,  rigid 
inconel  electrode  clamped  onto  a  tapered  nickel  plug,  vio-ch,  in  trim,  is 
brazed  to  the  tube  near  the  lower  end,  and  through  a  graphite  plug  to  a 
tantalum  electrode  at  the  vpper  end  of  the  test  section.  Tiie  centex’  of  the 
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Figure  16  4000”F  Test  Section  Assembly 
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tantalum  electiXKie  is  filled  with  graphite  powder  compressed  between  two 
graphite  retaining  rings. 

Pressure,  applied  to  the  upper  retaining  ring  by  means  of  a  tension^  » 

loaded  magnesium  piece,  insures  good  electrical  contact.  Since  the 
graphite  offers  little  resistance  to  slip,  the  tantalum  test  section  is  free 
to  undergo  thermal  expansion. 

Insulation  of  the  test  section  is  provided  hy  two  materials:  fully- 
stabilized  zirconium  oxide  powder  in  regions  of  high  ten5)erature  and  micro- 
quartz  fibrous  insulation  at  lower  temperature.  The  zirconium  oxide  powder 
is  packed  into  an  annulus  formed  by  a  l/2-inch  diameter  tantalum  tube,  which 
also  serves  as  a  radiation  shield,  and  a  2.5  inch  diameter  inconel  tube. 

Micix>-quartz  insulation  fills  the  remainder  of  the  pressxire  vessel,  Ihe 
insulation  reduces  the  radial  heat  loss  from  the  test  section  to  about  10 
per  cent  of  the  total  power  input  and  limits  the  maximum  teirg^erature  of  the 
pressure  vessel  walls. 

The  outside  wall  temperature  of  the  heated  test  section  is  measured 
at  six  axial  locations  by  a  radiation  pyrometer  sighting  through  the  view 
windows  provided  on  the  pressure  vessel.  Small  tantalum  tubes  pass  -through 
the  insulation  to  give  a  clear  view  of  the  test  section.  In  addition, 
wall  temperatures  are  also  measured  with  platinum-platinum- rhodim  themo- 
couples  welded  to  the  test  section.  These  thermocouples  are  limited  to 
temperatures  below  3100^F,  but  they  provide  a  valuable  check  on  the  pyrcsaeter 
readings  at  these  temperatures. 

The  helium  leaves  the  pressure  vessel  -through  a  water-jacketed  inconel 
tube  ^diich  forms  the  first  section  of  a  calorimeter.  The  calorLnater  cools 
the  exhaust  gas  to  a  temperature  that  can  be  measured  by  a  chromel-alumel 
"thermocouple  placed  in  a  mixing  chamber.  Additional  measuremen-ts  of  the 
gas  and  wa-ter  flow  rates,  and  of  -the.  entering  and  leaving  water  temperatures, 
allow  the  bulk  temperature  of  the  gas  leaving  the  test  section  to  be 
calculated.  A  second  cooling  section  is  ins-tailed  in  the  exhaust  line  as 

the  helium  is  still  too  hot  after  passing  through  the  calorimeter  to  flow  , 

immedia-bely  "through  the  exhaust  -valve. 
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